1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.00
## LDevsum 1.00 1.01
## dh0 1.10 1.26
## dh1 1.16 1.41
## dh2 1.05 1.13
## dh3 1.01 1.02
## dl0 1.04 1.06
## dl1 1.12 1.32
##
## Multivariate psrf
##
## 1.11
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1146.14620 | 22551.3845 |
| DIC3 | 1112.94318 | 22642.5526 |
| PWAIC | 23.35143 | 206.4733 |
| WAIC | 1124.71772 | 22656.2759 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.03
## LDevsum 1.00 1.00
## dh0 7.43 13.83
## dh1 16.90 30.91
## dh2 14.34 32.04
## dh3 8.50 16.56
## dl0 1.25 1.63
## dl1 1.10 1.27
##
## Multivariate psrf
##
## 20.1
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1387.28489 | 33774.328 |
| DIC3 | 1249.46624 | 33888.911 |
| PWAIC | 92.30564 | 167.363 |
| WAIC | 1319.99773 | 33893.412 |